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We introduce a novel class of field theories where energy always flows along timelike geodesies, 
mimicking in that respect dust, yet which possess non-zero pressure. This theory comprises two 
scalar fields, one of which is a Lagrange multiplier enforcing a constraint between the other's field 
value and derivative. We show that this system possesses no wave-like modes but retains a single 
dynamical degree of freedom. Thus, the sound speed is always identically zero on all backgrounds. In 
particular, cosmological perturbations reproduce the standard behaviour for hydrodynamics in the 
limit of vanishing sound speed. Using all these properties we propose a model unifying Dark Matter 
and Dark Energy in a single degree of freedom. In a certain limit this model exactly reproduces the 
evolution history of ACDM, while deviations away from the standard expansion history produce a 
potentially measurable difference in the evolution of structure. 



I. DUSTY FLUID WITH PRESSURE? 

How can one obtain dust from a scalar field? One can 
imagine a canonical scalar-field where the kinetic term is 
constrained to be equal to the potential. We can imple- 
ment this property by introducing a Lagrange multiplier, 
A, in the Lagrangian, 



We then find that the pressure is identically vanishing 
on all solutions and energy follows geodesies. This model 
describes the usual dust without vorticity. 

How can we obtain "dust with pressure" ? We can gen- 
eralise the above by adding some function of the scalar 
field and its derivatives to the Lagrangian, 

The constraint remains in effect and standard scalar-field 
dynamics are not restored. In fact, we will show that 
fluid elements in all such theories also always flow along 
geodesies, mimicking in that respect standard dust, yet 
the fluid has non- vanishing pressure. With this simple 
idea we have separated the notion that the pressure of 
the fluid is tied to the motion of a fluid element as is the 
situation in the usual case, e.g. radiation or cold dark 
matter. A parcel of such fluid will flow along geodesies, 
yet a manometer will record a pressure changing with 
time. 

In this paper, we introduce this new class of scalar-field 
models, which we will call X(p-fluids. These theories are 
described by an action containing two scalar fields, ip and 



*Electronic address: cugene.a.lim@gmail.com 
tElectronic address: ignacy.sawicki@nyu.edu 
■t Electronic address: alexander.vikman@nyu.edu 



A, where the latter plays the role of a Lagrange multiplier 
and enforces a constraint relating the value of the scalar 
field (fi to the norm of its derivative. This constraint 
forces the dynamics of the Xip-fluid to be driven by a 
system of two first-order ordinary differential equations, 
one for the fleld ip, the other for the Lagrange multiplier. 
As a consequence, there are no propagating wave-like de- 
grees of freedom and the sound speed for perturbations 
is exactly zero irrespective of the background solution. 
However, the initial-value problem still requires the spec- 
iflcation of two functions on the initial time slice. Thus, 
effectively, a single dynamical degree of freedom remains. 

Provided that the derivatives of the scalar field (p are 
time-like, the system can be interpreted as a perfect fluid. 
However, for a Xtp-fluid given by a particular action, the 
relation between the pressure p and the energy density 
£ is solution dependent. We show that an arbitrary ef- 
fective equation of state, including phantom ones, can be 
obtained by choosing the form of the Lagrangian appro- 
priately. In addition, for all Xip-fiuids, there always ex- 
ists a region in their phase spaces in which the Xip-fluid 
is effectively pressureless. We will exploit this feature to 
model the evolution of the cosmological background from 
matter domination through to the acceleration era as be- 
ing driven by the dynamics of a single degree of freedom 
provided by the Xip-fluid. 

The key novel aspect of this class of theories is that 
the four-acceleration is always zero, even when the pres- 
sure does not vanish. This is a result of the constraint's 
eliminating those fluid configurations where the pressure 
has a gradient orthogonal to the fluid velocity. 

Motivated by these properties, we will use the Xip-fluid 
to frame the problem of the dark sector in cosmology in 
a unified manner. The existence of the dark sector in 
the Universe's energy budget is now established beyond 
reasonable doubt. The standard model of cosmology, 
ACDM, splits it into two constituents: cold dark mat- 
ter ("CDM") — a pressureless fluid ("dust") which clus- 
ters allowing baryonic structures to form in its poten- 
tial wells and detectable to this day in the form of ha- 
los around galaxies — and dark energy ( "DE" ) — a form of 
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energy that appears to be smooth and to have an equa- 
tion of state close to a cosmological constant. This di- 
chotomy of phenomenology has made it difhcult to build 
a compelling model which would treat the two dark com- 
ponents in a unified fashion. Nonetheless, some models 
exist in the literature: [1-13]. It is not hard to see that, 
given a fiuid which clusters like dust yet has arbitrary 
pressure, we can construct such a unified model — which 
we call Dusty Dark Energy ("DDE"). We present the 
main results of our application of Xip-fluid to such a cos- 
mology beneath. We refer the reader to the full analysis 
in the main body of the paper, section V. 



A. Summary of Cosmological Results 

We have studied cosmological perturbations in the case 
when an arbitrary Xip-fluid dominates the Universe. We 
have derived the closed-form equation for the evolution 
of the Newtonian potential $ which turns out to re- 
cover the standard result for general hydrodynamics in 
the limit of vanishing sound speed. This evolution is 
determined by background expansion history only and 
in the limit of the ACDM expansion history the evolu- 
tion of perturbations is exactly as in the standard case. 
We have also written down an action for perturbations 
which explicitly shows that there are no ghosts in this 
theory when the equation of state for the Xip-fluid is non- 
phantom. We demonstrate that on a classical level our 
model can cross the phantom divide without singular- 
ities while linear perturbations continue to evolve sta- 
bly; however, the perturbations do become ghosts at this 
point, hence making the system unstable (in particular 
quantum-mechanically) when interactions are taken into 
account. In this paper we put the investigation of that 
instability and issues related to the possible strong cou- 
pling scales aside. 

We consider a universe comprising only radiation and 
the Xip-fluid which will describe both CDM and DE. To 
illustrate more concretely the phenomenology we have 
focused on a specific family of models which is parame- 
terised by wnn — the equation of state of the Xip-fluid in 
the asymptotic future. Given that the initial values for 
the Xip-fluid are chosen appropriately, the radiation be- 
comes subdominant while the Xip-fluid is still far off its 
final attractor (given by wnn) and evolves approximately 
like dust, giving an epoch of matter domination. The 
duration of this epoch is also determined by the initial 
values. 

In the limit — > — 1 this family of models recovers 
exactly the background evolution and growth of structure 
of ACDM. However, if Wfin ^ —1, the evolution of the 
background and perturbations differs from a "wCDM" 
model comprising cold dark matter and dark energy with 
a constant equation of state. We illustrate the evolution 
of the effective equation of state for the dark sector in 
a selection of different Wfin-cosmologies in Fig. 1. We 
find that the transition from matter-domination to dark 
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FIG. 1: Time evolution of the total effective equation of state 
for the dark sector. The black solid line represents the evo- 
lution in ACDM which is identical to that of the lUfin = — 1 
model. Models with final equations of state 1 -I- Wfln > begin 
to deviate from matter domination earlier and the transition 
is slower than ACDM. The opposite is true in phantom mod- 
els. The evolution is normalised such that the equation of 
state at a = 1 matches the best-fit result for the ACDM cos- 
mology as determined by WMAP7 results, wo = —0.74 [14]. 



energy domination differs from that of wCDM (Fig. 2). 
We compare the growth of linear perturbations with that 
of ACDM in Fig. 3. We find that models with 1 + Wfin > 
exhibit a growth factor suppressed by a few tens of 
percent. Also the Newtonian potential here is lower by 
a few percent than in ACDM, which would decrease the 
integrated Sachs- Wolfe effect. 



B. Open Questions and Future Directions 

Finally, we mention open issues in this setup which 
remain to be addressed 

• The nature of caustics: It is well known that in non- 
canonical field theories caustic can develop, e. g. 
see Ref. [15] on caustics in Sen's string-theoretical 
tachyon matter [16, 17] and Ref. [18] on caustics 
in the ghost condensate [19] and the discussion in 
Refs [20, 21] on caustics in Hofava gravity [22]. We 
expect that the Xip-fluid will develop caustics. The 
question is how to interpret the multivalued regions 
once this occurs. 

• Can Xip-fluids virialise? If the Xip-fluid is to really 
model non-linear structure, it must be able to form 
static and stable configurations (e.g. halos). 

• What is the origin for the initial conditions for 
Dusty Dark Energy introduced in section V? Could 
a more generalised setup provide a solution for the 
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FIG. 2: Comparison of the derivative of the effective equa- 
tion of state for Dusty Dark Energy (Eq. (41)) with that for 
a dark matter plus dark energy with a constant equation of 
state, wCDM, (Eq. (43)). The magnitude of the derivative de- 
termines the duration of the transition between matter dom- 
ination and the acceleration era. The left panel shows that 
for Wfin > — 1, the transition in the DDE model is more rapid 
than the corresponding wCDM model. On the other hand, 
for phantom Wfin this transition is slower than for the corre- 
sponding ii;CDM model, as shown in the panel on the right. 
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FIG. 3: The comparison of the total growth of perturbation 
amphtude between Dusty Dark Energy and ACDM. The evo- 
lution of the Newtonian potential, $, is determined by Eq. 
(25) and deviates by a few percent from its ACDM values by 
a few percent. This evolution will affect the strength of the 
ISW signal in the CMB. On the other hand, the evolution of 
the density perturbation on subhorizon scales is determined 
by Eq. (22) and is affected much more strongly. 



coincidence problem? In our model, we require that 
at some point during the radiation epoch the en- 
ergy density of the XLp-fluid be equal to the energy 
density of CDM: we have no alternative to the stan- 
dard dark-matter freeze out scenario which would 
produce this in a natural manner. 

• What is the Hamiltonian structure of this theory? 
How to quantise it? See, for example [23]. Is the 
structure of this theory stable as a result of ra- 
diative corrections: would a kinetic term for A be 
generated? 

• What is the strong-coupling scale for perturba- 
tions? We should note that our model is rather 
similar to a potentially singular limit of Hoi^ava 
gravity [22]. There it was found [21] that, contrary 
to our case, the sound speed for cosmological per- 
turbations becomes imaginary, and that the strong 
coupling scale may be extremely low. 

• What is the rate of instability for the phantom 
case? Is this instability catastrophic as it is usu- 
ally for ghost degrees of freedom? It is possible 
that the absence of propagating wave-like degrees 
of freedom may change the standard picture [24- 
28]. 

• The action formulation of this theory allows us to 
consider couplings to standard model fields. It is 
natural, for example, to consider Lorentz violation 
in this framework in analogy to Einstein aether the- 
ories [29-31]. 

• Can this theory be a low-energy limit of a more 
fundamental theory? 

• Can our conceit with the constraint be usefully 
extended beyond scalar fields to theories with 
fermions, vector fields or many degrees of freedom? 



II. DYNAMICS IN GENERAL CURVED 
SPACE- TIME 

Let us consider a scalar field theory given by the action 

1 



5 = 



where the field A is a "Lagrange multiplier" and does not 
have a kinetic term, while 



1 



is a standard kinetic term for the field K {(p, X) is an 
arbitrary function of X and (/?, while n (tp) is an arbi- 
trary function of the scalar field </?. In hydrodynamical 
language, ip is one of the velocity potentials and does 
not necessarily have to represent a fundamental degree 
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of freedom; plays the role of the specific inertial 

mass [32]. Let us further assume a standard minimal 
coupling with gravity. 

The equations of motion are 



1 SS 



-A/i^^- V„ (AV"(/j) =0. 



(1) 
(2) 



Here and throughout the paper we denote partial deriva- 
tives by subscripts. The Lagrange multiplier imposes a 
non-holonomic constraint^ on the dynamics of the fields, 
hence the action cannot be written purely in terms of Lp. 
In the case of time-like derivatives: X > 0, similarly to 
k- essence [33-35], we can introduce an effective 4- velocity 



(3) 



where in the last equality we have used the "constraint" 
Eq. (1). Further, we can see that the corresponding ef- 
fective four-acceleration always vanishes. 



(4) 



Here and throughout the paper we use the notation ( ) = 
u^Va ( ) for the derivative along u". Thus, on equations 
of motion, are tangents to time-like geodesies. It is 
also convenient to write Eq. (1) in this form 



The general solution of Eq. (5) is 

^(a;") = /(r-^(x)) , 



(5) 



(6) 



where r is a time parameterising the congruence of 
geodesies given by u", / is a general function solving 
Eq. (5) and (x) is an arbitrary function of spatial co- 
ordinates X in the hypersurface normal to u". One can 
consider ip to be an intrinsic clock for our system, since 
Eq. (6) defines a time reparameterization. 
The energy-momentum tensor (EMT) is 



Tap — 



5S 



-g Sg 



ap 



= {Kx + A) Va^Vpif - Kgap , 



which is of the perfect-fluid form. Using hydrodynamical 
notation the EMT can be rewritten as 



P5a/3 : 



Tap = {e+p) UaUp 

with the energy density given by 



(7) 



^ As a technical point, the system possesses two second class con- 
straints, one which is primary and the other secondary. 



where into Kx and K we substitute the constraint X = 
^IJ,^. The pressure is a function of the intrinsic clock ip 
only 



p{ip) = K [ip 



(8) 



This is the key feature responsible for the absence of ac- 
celeration, Eq. (4), since the gradient of the pressure is 
always parallel to . This can be explicitly seen by con- 
sidering the conservation of the EMT: 

V„T"^ = {e+p + 0ie + p))u^ + {e+p) -V^p = 
= {e + e{e+p))uf^ , 

where we have used the fact that pressure does not have 
gradients orthogonal to and = VaU"^ is the expan- 
sion of a congruence of geodesies under consideration. 
For perfect fluids in the usual case, the divergence of the 
EMT has two components: one parallel to the velocity 
describing the conservation of energy and one parallel 
to the four- acceleration . Here the latter is identically 
zero and so the conservation of the EMT reduces to the 
conservation of energy. 



e + e{e+p) = Q. 



(9) 



Then by choosing K and /i appropriately one can arrange 
for a general evolving equation of state 



wx 



P 

e 



Despite this, the energy flux of the Xip-fluid, T""up = 
eu", always follows time- like geodesies, as is the case for 
perfect fluids in the absence of pressure. Note that in 
FRW universes, the energy of an arbitrary fluid flows 
along time-like geodesies as well, provided that the fluid 
configuration be bound to respect the homogeneity and 
isotropy of the FRW spacetime. This is of course not true 
of inhomogeneous perturbations in the standard case. 

Now, let us refocus on the equation of motion. Using 
the constraint Eq. (1) we can re-express the gradient of 
the kinetic term and the expansion in terms of ip, 



VaX = ^J.p.^VaP', 

9 — \7aU" = ^"^□(/9 — /Xy . 



(10) 



The system of equations of motion Eq. (2) and Eq. (1) 
can be written as 



A = -^"^ {e^pi + (e + p) ( 



(11) 
(12) 



where for the partial derivative of e^p we differentiate Eq. 
(7) to obtain 

= ^l^lp {fi^Kxx +Kx + 2A) + fi^Kx^p - . (13) 
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The equations of motion, Eqs (11) and (12), for our sys- 
tem have been reduced in this way to two first-order or- 
dinary differential equations. They then have to be sup- 
plemented by the Landau-Raychaudhuri equation 



e 



1 



a/3 



for the expansion 9 and similar equations for the shear 
(Tq,/3 or, equivalently, by Einstein's equations'^. 

From the above system of equations of motion and the 
4-velocity Eq. (3) it follows that the Cauchy problem lo- 
cally has a unique solution depending on two functions 
if (x) and A (x) given on an initial spacelike hypersurface 
S. Moreover, from Eq. (11) and Eq. (12) it follows that 
the solutions propagate along time-like geodesies. There- 
fore in the future of the geodesic 7^ starting at Xj on E 
the solution only depends on the initial values of Lp (x^) 
and A (x^) at the point Xj — the sound speed Cg is identi- 
cally equal to zero for all configurations of Xip-fluid. This 
means that there are no wave-like dynamical degrees of 
freedom. Adding dynamical gravity does not change the 
picture. 

As we have already mentioned in the introduction, the 
global Cauchy problem may be ill-defined for some initial 
data since caustics may develop. However, this problem 
is not unusual for non-canonical field theories. 

Finally we note that, since £\ = fj? ^ 0, we could use 
the pair (</?, e) instead of (ip, A) as independent variables. 
In that case, one has to use energy conservation Eq. (9) 
instead of Eq. (12) as the second equation of motion. 



III. DYNAMICS IN COSMOLOGY 



A. Cosmological Background 



where Mpi = (SttGn) is the reduced Planck mass. 
The phase space for this system of two first-order ordi- 
nary differential equations is (A,(y5). 

For completeness we also present the second Fried- 
mann equation 



H = - 



A*' {Kx + A) 



2M2j 



2Af2j 



(16) 



The presence of an external energy density of some 
other cosmological fluid brings with it the correspond- 
ing change to the Friedmann equations, Eqs (15), (16) 
but no changes to the equations of motion apart from 
those implied by the change in H. 



B. Cosmological Perturbations 

Symmetry considerations imply that our scalar field at 
linear order only contributes to the scalar part of the cos- 
mological perturbations. In Newtonian gauge, the metric 
for scalar perturbations is 

ds^ = (1 + 2$) - (t) (1 - 2^-) dx^ . 

The absence of anisotropic stress, <& = vji^ simplifies the 
perturbed Einstein equations to (see e.g. [36]): 



fetot 



$ + 4i/$ + {2H + 3H^) $ 



Sptot 
2A/2^ 



(17) 
(18) 
(19) 



In the case of a background cosmology of a pure Xtp- 
fluid in the FRW universe with the metric 

ds^ = dt^ - (t) dx^ , 

our equations of motion take the form 

v'=^^iv), (14) 

X = - fi-^ {s^fi + 3H{e+ p)) , 

where we have used 6 — 3H = 3a/a. Alternatively, we 
can use energy conservation Eq. (9) instead of the equa- 
tion for A. To close this system we use the Friedmann 
equation-^ 



3M2j 



31/2; 



(15) 



^ Note that the rotation 0;^^ is zero because of Eq. (3). 

^ Here, for simpHcity, we have assumed a spatially flat case. 



where A — didi and the (5utot||j: is the potential (scalar) 
part of the i component of the total four-velocity. Per- 
turbing the equations of motion gives: 



for the constraint Eq. (1) and 



(20) 



5e- 



A f Sip 
3$ + ^ ' ^ 



{e+p) + '3H{Se + 5p) = 0, (21) 



for the perturbations of the energy density. See Appendix 
A for the derivation. This equation is standard for hy- 
drodynamical matter [36, p. 312, Eq. (7.105)]. The non- 
standard input of our model is that Sp = p^Sip and that 
Eq. (20) describes the evolution of the velocity poten- 
tial V — fi~^Sip for time-like geodesies. Note that in the 
above equations, we have not assumed the domination 
of the Xtp-fluid. Thus, equations (20) and (21) can be 
used to follow the dynamics of the linear perturbations 
of Xip-fluid through the entire history of the universe. 
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For future discussion it is helpful to re-express the above 
as an equation for the relative perturbation = Seje: 



4 = - 'SH^Sifi^ 



(22) 



3$ 



A (5lp 



a 



[l + wx] 



Combining the above with the perturbed constraint Eq. 
(20) and obtaining the Newtonian potential $ through 
the perturbed Einstein equations (17), (18), (19) closes 
the system and allows us to describe the evolution of the 
perturbations in the Xtp-fluid in general case. 

When the Xip-fluid dominates the energy density and 
the perturbations, a significant simplification occurs since 

^Utot||i =H^^diSip, 
Sstot =£ipSip + f/6X , 

Using these expressions for perturbations, we can first 
integrate Eq. (18) to obtain 



Sip — 



2M2j^ 



P 



$ + = - 



H 



where in the last equality we have used the second Fried- 
mann equation (16). Combining this with Eq. (19), we 
can obtain a closed expression for $ 



H 



(23) 



H^<^ = 0. 



Thus, the Newtonian potential always evolves as 

<i>(<,x)=/i(i)Ci(x)+/2(i)C2(x) , 



(24) 



where Ci (x), C2 (x) are arbitrary functions of the spa- 
tial coordinates while /i {t) and /2 (t) are solutions of 
the homogeneous ordinary differential equation Eq. (23). 
Using the solution Eq. (24) we can find all other quan- 
tities. In particular, substituting this solution for <I> into 
the Poisson equation (17), we obtain Se. The separation 
of variables in the solution (24) implies that, as we have 
discussed above, the sound speed is zero. 

It is convenient to introduce dimensionless time (e- 
folds number) iV = In a and rewrite the differential equa- 
tion for the Newtonian potential $ in terms of N 



+ $' 4 



H' 



3 + 2^+r)$ = 0, (25) 



where we have introduced a dimensionless correction to 
the equation for standard pressureless dust arising from 
the perturbation of pressure in our model. 



2M2jiJ'i/2 



(26) 



with ( )' = On ( ) the derivative with respect to the num- 
ber of e-folds. Using the Friedmann equations, F can be 
written in geometrical terms as 



F = -- 



El 



-3. 



This particular combination vanishes for a background 
expansion history mimicking that of ACDM. Therefore, 
in that limit, linear perturbations in our model will be- 
have identically to dust in the presence of a cosmological 
constant. 

Further, similarly to the standard case one can intro- 
duce a new variable 



Q 



-H 



such that Eq. (25) takes the form of an oscillator with 
time-dependent frequency 



e 



Q = , where 6 



H 



(27) 



Our variable Q is a redefinition Q cx uyaH of the stan- 
dard variable 



(28) 



given in [36, p. 302, Eq. (7.63)]. This redefinition is 
caused by our choice of the time coordinate. The equa- 
tion of motion for u is 



dlu 




u = 0, 



where r) is conformal time [adrj = dt) and 




(29) 



(30) 



It is easy to check that equations (27) and (29) are equiv- 
alent. Thus, cosmological perturbations of the Xip-fluid 
reproduce the standard result for hydrodynamical matter 
in the limit Cg = 0. We could have guessed the equation 
of motion (29) and the variable u from the very begin- 
ning because neither the final expression for u nor the 
formula for 9 explicitly involve the sound speed. How- 
ever, note that the derivation of these results presented in 
[36, p. 302] uses in an essential way the standard hydro- 
dynamical formula for the sound speed = p/i which 
is clearly absolutely inapplicable for the Xip-fluid. Curi- 
ously, u and 9 are given by the same formulae (29) and 
(30) in the case of k- essence where the derivation uses a 
sound speed given by Cg — px/^x 7^ p/e (see [36, 37]). 

From equation (27), one can see that one of its solu- 
tions is Q cx 0, which translates to one of the modes for 
the Newtonian potential 



$ cx 



H 
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Taking the derivative gives the "instantaneous power 
law": 



In 



H 



'iwx \ 



2 / 



demonstrating that for wx > —5/3 this solution repre- 
sents the decaying mode. Following the discussion in [36, 
p. 303], we find the second mode of the solution from the 
Wronskian. Integrating this result by parts allows us to 
write the full solution as 

$ = Ci (x) + -C2 (x) ~ Ci (x) - (31) 
a a J H 

This solution is valid on all scales, whereas in the stan- 
dard case when Cg 7^ it is only applicable on superhori- 
zon modes. 



1. Perturbations Around Scaling Solutions 

Let us analyse this general solution further in the case 
of a cosmology dominated by a Xip- fluid with a constant 
equation-of-state parameter wx 7^—1 (we demonstrate 
how to construct one such class of models in section 
IV C). In such a case, the final term in the expression 
above is constant, and our solution for wx 7^ —5/3 is 

$ = C'i(x) + C2(x)a-(5+3"'^)/2^ (32) 

while for the special case of wx = — 5/3 the solution is 

*-5/3 = C'i (x) + C'2(x) In a. 

Thus, for all constant wx > —5/3 the Newtonian poten- 
tial is constant up to a decaying mode. Substituting this 
solution (32) into the Poisson equation (17), we obtain 
for the density fluctuations in the cases wx 7^—1 and 
Wx ^ -5/3: 

2 A 

5- = -2$ -I- - 2$' = 

3 [aHf 

= -2Ci +3C2(l + u'x)a"<'+^"'^^/' + 
For subhorizon modes, k ^ aH , we have 



l+3u)x 



for Wx > 



while for superhorizon modes, k <^ aH , we have 



3 

e„»aH-a'''"''-''^' for .,^<-^ 



i^e)k^aH const 



for Wx > 



('5e),«,^'-a-(^+3--)/2 for wx< 



It is illustrative to compare these results with those for 
standard cosmological fluids. In particular, for the ultra- 
relativistic equation of state, wx — 1/3, the Newtonian 
potential for the Xtp-fluid remains constant up to a de- 
caying mode at all scales, while subhorizon density per- 
turbations grow as Se ^ a^. This should be contrasted 
with standard radiation which causes both the potential 
and the density perturbation to decay and oscillate once 
the mode becomes subhorizon. 



2. Phantom Behaviour and w — —1 Crossing 

As the general solution (31) implies, the perturbations 
do not have catastrophic instabilities even for phantom 
[38] equations of state, with wx < —1. The \ip- fluid 
framework allows one to realise such scenarios easily. 
However, Eq. (26) implies that F has a singularity when 
H' = 0, i.e. when wx — —1- Can a \ip- fluid evolve 
through this singularity, from a standard to a phantom 
equation of state? 

The answer is provided by the analysis of Eq. (25) in 
the vicinity of the singularity. The derivatives of $ can 
only remain finite, if the singular terms in the equation 
cancel, i.e. if $' -|- $ = at the phantom-divide-crossing 
point. Indeed, it can easily be checked that the general 
solution (31) always satisfies 



<f'+$=^(c2(x)-ci (x) 



so that automatically $' -|- $ = when H' = 0. There- 
fore, we have also shown that in our model there are 
no classical catastrophic instabilities associated with the 
crossing of the phantom divide. 

Note that the Xip-fluid contains only one degree of free- 
dom but the action cannot be written exclusively in terms 
of this degree of freedom in a generally covariant local 
form. Therefore this possibility of smooth crossing of 
the Wx = — 1 barrier does not contradict to the state- 
ment proved in [39] and rederived in different ways later 
in [40-46]. The Xtp-fluid provides a working example 
of the so-called Quintom scenario from [47], see also re- 
views [48-50]. Further, abandoning another assumption 
from [39] , that coupling to gravity is minimal, allows one 
to have a classically stable crossing of the phantom di- 
vide in scalar-tensor theories [51-55]. For other single- 
field options see [56-58]. For more on phantoms see e.g. 
[8, 59, 60] and references therein. 

One can obtain the equation of motion (29) from the 
action 




(33) 



Both quantities u and 9 (or Q and &) are defined up 
to a constant factor. In particular, this factor can be a 
complex number e. g. the imaginary unit i. Note that 
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the definitions (28), (30) which we have used imply that 
u and 9 are real, provided that the Null Energy Condi- 
tion for the background holds i.e. H < 0. As usual, 
the sign of the action (33) is such that for H < {) \i 
has a positive definite kinetic term. When H — 0, both 
quantities u and diverge. However, as we have already 
demonstrated, the evolution of the physical quantity $ 
is smooth through H ^ 0. This also means that the cur- 
vature invariants are smooth through wx = ^1 crossing. 
Thus after the crossing u and 6 become pure imaginary 
or in terms of real fields, the action changes the overall 
sign and the kinetic term becomes negative definite. 



IV. EXAMPLE Ai^-FLUID COSMOLOGIES 

A. A</p-Dust with Cosmological Constant 

Let us consider the simplest case: — and fi — 
const. In that case, we have 

p {fi) = K and e = [Kx + X) - K , 

where K (/i) = const and Kx (a*) = const. For the equa- 
tions of motion, Eq. (14), we have 

X = -3H{Kx + X) , 



with the solution 



(Kx + X) 



where Eq is a constant of integration. Thus the system 
behaves as 



p — const and e — Equ 



so that the energy-momentum corresponds to a mixture 
of a cosmological constant A = —K (/i) of either sign 
and pressureless dust with the energy density eo today. 
Since the background is that of ACDM cosmology, by the 
argument of section IIIB, the evolution of $ and 6^ also 
exactly reproduces the standard results. 



B. A(/3-Dust 

If we take K — and an arbitrary /i (ip) we obtain 
a A(/3-fiuid which mimics pure dust, p — 0, with energy 
density e = pi^X. Indeed from Eq. (14) we have 



Xn^ 



- — 1[iiiX 



iHe, 



with the standard solution e = e^a^^ . Note that there 
is a degeneracy in ip {t) . In this setup for different /i {ip) , 
the same evolution of the energy density e [t) corresponds 
to different ip [t). Again, by our discussion of cosmologi- 
cal perturbations in section IIIB, the evolution of $ and 
exactly reproduces the results for the standard dust- 
dominated universe. 



C. Xp-Fluid Possessing a Scaling Solution 

Moving beyond a constant /z, let us consider a class of 
models with 



K = aX , where 
/I = Mocxp (-— ) , 

where the mass scale for p is 



cr = ±1 



^/q^Wfin 

3 1+ Wfin ' 



PI • 



(34) 
(35) 



(36) 



In the following, we will show that, in this class of mod- 
els, the dynamics of the A</3-/Zwrf-dominated cosmological 
background has a fixed point with a constant equation of 
state Wfin, where Wf^■c^ can have either sign and can even 
be phantom-like.*. Further, we will show that this fixed 
point solution is an attractor provided Wfjn < 1. 
From Eqs (8), (7) we obtain 



p - 



-ji and e = ^ V 2 ^ "^7 ' 



so that the instantaneous equation of state, 

1 

Wx = 



1 + 2ctA 



(37) 



(38) 



is determined by the value of A. From this equation it 
follows that, if wx = const, then 



X 



— 1 ) = const . 



(39) 



In particular, for values of wx corresponding to an accel- 
erating expansion, A^,^ is a number of order one. Mean- 
while, the exponential form of ^ {p), Eq. (35) implies that 
the evolution of /i has a very simple form. 



TO 

n = — and 
^ t 



TO 



(40) 



where we have chosen constants of integration in such a 
way that the pressure p of the Xp- fluid is singular exactly 
at the Big Bang, t — Q. 

We can now rewrite the equation of motion Eq. (14) 
combined with the Friedmann equation for a universe 
containing solely the Xp-fluid, to obtain the equation of 
motion for the equation of state 



w 



X 



iwx 1 + Wx 



Wx 
Wfin 



^1 + Wf^,,] 



(41) 



It is easy to see that Wfin is the fixed point of this equa- 
tion, and therefore also of the equation of motion for A. 
In the limit of Wfjn — > —1, this equation reduces to the 



- Obviously we cannot realise uign = in this setup. 
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evolution of w for ACDM. Therefore, as wnn approaches 
that hmit, the transition from matter domination to dark 
energy domination becomes indistinguishable from that 
in ACDM. 

Let us now consider the stability of this fixed point. 
Linearising the above around wx — wnn, we obtain 

3 

^'^'x = 2('^fin ~ l)Swx ■ 

From this result, it follows that in an expanding universe, 
for uifin < 1) the equation of state approaches wan and 
A approaches the fixed point Au,fi„. We call this solution 
the w-attractor. 

Eq. (41) can actually be solved, albeit implicitly, al- 
lowing us to obtain the scale factor a as a function of the 
equation of state parameter wx- 



3(u.fi„-l) 



y^Wfin/wx - lUfin 



(42) 



therefore resembles a matter-dominated epoch, eventu- 
ally approaching an era with a constant wx = wnn- 

Approximating Eq. (41) for small wx we can calculate 
how the equation of state evolves during this period: 



3wx 



Wx 



A 



-,-3 



Turning to perturbations during this matter-domination 
era, the pressure corrections Eq. (26) can be written as 



r (wx) = 3wx 



Wx f 1+ Wan 



WRn V 1 + Wx 



-9/2 



(44) 



where we can explicitly see that F -> as wnn — J> — 1. For 
parameter values motivated by dark energy, 1 1 -I- w^n \ ^ 
1, these pressure corrections to the evolution of the po- 
tential $ are subleading until wx approaches its attractor 
value, since in Eq. (25) 



wx+T)+ {-3wx + F) $ = . (45) 



where oq is a constant of integration. 

It is instructive to compare the evolution of the Xtp- 
fluid given by Eq. (42) or Eq. (41) with the case of a 
mixture of dust and DE with a constant equation of state 
w = Wfin- For iffin < the evolution of the universe has 
the same late-time asymptotic as in the Xip-fluid case. 
From the Friedmann equations (15) and (16) we obtain 



^dark = 3Wdark (wdark " Wfin) , 



(43) 



instead of Eq. (41), where we have denoted the total ef- 
fective equation of state as Wdark- These two equations 
(41) and (43) only coincide in the limit of w^n — 1- 
For illustrative purposes we also present the solution for 
(43) in the form similar to Eq. (42): 



Wdark 



- 1. 



From this analysis then, the evolution history of the Xip- 
fluid cannot be reduced to the evolution of a mixture of 
DE with a constant equation of state parameter wgn and 
dust, providing a potential observational probe for such 
theories. 

In the following, we present cosmologies in the two lim- 
its where the Xip-fluid dominates the total energy density 
of the universe and when it is subdominant to some other 
matter fluid. 



1. Dominant Xif -Fluid in FRW Universe 

Let us again consider the model given by Eqs (37) and 
(34), however, this time in the case with A ^ 1, far off 
the w-attractor. Here, the system is effectively pressure- 
less, by virtue of Eq. (38). The energy density scales 
as since wx ~ 0. The evolution off the attractor 



So, while off the attractor, there is no significant devia- 
tion in the growth of structure from standard considera- 
tions. 

On the attractor, wx = Wfin, we obtain 

F = 3Wfin , 

and Eq. (45) confirms what we found in section III B for 
the general model: for w^n > —5/3, $ = const, while 



,l + 3u)tin 



is the growing mode inside the horizon. 



2. Subdominant Xtp-Fluid in FRW Universe 

Motivated by the discussion above let us still consider 
theories with 



= ^0 exp 



(r-) 

\ m/ 



where m — const and may, for example, be given by Eq. 
(36). Using the result of Eq. (40) we obtain the equation 
of motion for A, 



A = i (ct + 2A) 



3H{(7 + X) 



(46) 



In the history of the universe there were at least two 
stages with the background equation of state, Wb, ap- 
proximately constant, namely the radiation-dominated 
epoch — with Wh — ^ — and the matter-dominated era — 
with Wh — 0. Let us then consider the dynamics of Eq. 
(46) when the matter content of universe consists mostly 
of some fluid with a constant equation of state parameter, 
Wh = const. In that case, we have 



2/3(«.b + l) 



and consequently H = 



1 



3(wb + 1) t 
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here would correspond to the age of the universe today, 
when a = 1, if Wb = const. On this background, Eq. (46) 
for the subdominant Xip-fluid can be expressed as 



dA Wb - 1 , 2wb , 
a H -A . 



d In t Wb + 1 Wb + 1 
This equation has the following solution 



(47) 



A (t) =A^, + = A,„, + Ca^'"^ , for ^ , 

t 



Xit)^- In 



for W]^ — , 



where X^, is the fixed point solution given by the formula 
(39) and C and ti are constants of integration. We can 
therefore see that, for Wh 7^ 0, the Xip-fluid corresponds 
to a mixture of a fluid with the same equation of state 
Wh as the background and, in addition, dust 



P =77"^ t = 



a m 



2 V^o 

2 



-3(l+u;b) 



while for a dust-like background, w\y = 0, the Xip-fluid^s 
hydrodynamics obeys 



In the late-time asymptotic, the equation of state of the 
X(p-fluid approaches that of the background provided 
that Wb < 0. In particular, this tracking behaviour means 
that in an inflationary background when Wb ~ — 1, the 
Xip-fluid does not redshift away but instead survives in 
form of an effective cosmological constant. 



V. DUSTY DARK ENERGY 

In this section, we will present a unified dark matter 
and dark energy model using a single dynamical degree 
of freedom, which we call Dusty Dark Energy. 

The discussion of the constant-w model presented in 
section IV C is suggestive: a model with /i((y5) of the 
form given in Eq. (35) has a background evolution that 
evolves as dust when it is far off its attractor and even- 
tually settles on a constant equation of state determined 
by the value of the mass-scale parameter m. This picture 
is not altered when the Xip-fluid is subdominant during 
a radiation-domination epoch: the dust-like component 
will eventually overwhelm the redshifting radiation re- 
sulting in a period of expansion similar to matter domina- 
tion with its duration determined by the value of A when 
radiation becomes subdominant. Once A approaches a 



number of order unity, the transition to dark-energy dom- 
ination will occur at a rate similar to the usual transition 
in ACDM and eventually the expansion will settle on a 
constant (and arbitrary) equation of state. This back- 
ground evolution provides an expansion history that is 
close to ACDM and yet does not explicitly contain sepa- 
rate dark-matter and dark-energy components: only ra- 
diation and our Xip-fluid are necessary. 

Beyond describing the background dynamics, we have 
to show that the perturbations in the fluid evolve simi- 
larly to those in ACDM. This is necessary not only for 
the growth of large-scale structure during matter domi- 
nation, but also during the radiation-domination epoch: 
the anisotropies in the cosmic microwave background de- 
pend on the existence of appropriate gravitational poten- 
tials driven by dark matter. 

In section IV C, we have shown that during matter 
domination and at late times perturbations evolve simi- 
larly to the CDM perturbations in ACDM. The remain- 
ing piece is to prove that the perturbations evolve in a 
dust-like manner even when the energy density is domi- 
nated by radiation, which we will show as follows. 

During domination by radiation — or any external fluid 
in general — the Newtonian potential, is driven by per- 
turbations in that fluid. The subdominant components 
then respond to this potential and evolve according to 
equations of motion arising from the conservation of their 
EMT. For dust, the standard result, written in a sugges- 
tive manner, is 



Sip, 



dust 



(48) 



^dust 3$ 



A / Sipdust 
2 1, 



a 



In the case of a subdominant Dusty Dark Energy, the 
equivalent equations, Eq. (20) and (22) tell us that the 
eventual evolution to a matter-dominated era requires 
that A ^ 1, implying wx ^ 1/2 A. Therefore we can 
approximate 



Sip 



A 

3$-t- ^ 



-3Hm-^X-^S(p- -HX- 
^ 2 



(49) 



(50) 



All we now need to do is to show that the terms additional 
to the equations for CDM are negligible. For A ^ 1, 

4 - H5^ > HX^^Se ■ 

From Eq. (49) we can obtain the estimate 

/i<i> 



Sip 



H 



Then, the pressure correction in Eq. (50) is small com- 
pared to the Newtonian potential providing 
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or, in other words, if 



On the other hand 



1 > 



mHX 



(51) 



11 (1 + Wfin) \/-px 



mH 



PI 



(1 + Wfin) 



-Px 

Piad 



Since we assume that the DDE is subdoniinant, = 
3prad > \px\ and, by construction, |l + Wfin| ^ 1- 
Therefore all the additional terms in the evolution equa- 
tions are negligible, and the perturbations in the DDE 
obey the same equations for evolution as standard dust 
during radiation domination, demonstrating that the 
DDE model provides for a viable cosmology with our 
fluid playing the role of both cold dark matter and dark 
energy. 

The duration of the matter-domination era is deter- 
mined by the initial values of A(io) ^ 1 for any given 
initial ip{to). The former determines the time when the 
DDE turns over from dark-matter-like to dark-energy- 
like, while the latter sets the time where it begins to 
dominate over radiation. Since we have to tune both val- 
ues, we have not provided a solution to the coincidence 
problem. On the other hand, this model is a "minimal- 
ist" description — as we require only two initial conditions 
for set the transitions of the two epochs. 



VI. SUMMARY 

In this paper, we introduced a novel class of field theo- 
ries with a single dynamical degree of freedom which have 
a perfect-fluid interpretation. The key feature of this the- 
ory is that its fluid velocity flows along geodesies — hence 
mimicking "dust" in this respect. On the other hand, un- 
like a standard cold-dark-matter fluid, it carries pressure 
parallel to its fluid velocity. In cosmology, this pressure 
affects the expansion history. 

This sleight-of-hand is achieved by means of a "La- 
grange multiplier" field, employed to constrain by equa- 
tion of motion the above-mentioned behaviour of the fiuid 
velocity vector. Our system then consists of two first or- 
der equations of motion, and hence effectively a single 
degree of freedom. This dynamic cannot be reproduced 
by usual scalar field theories such as k- essence or higher 
derivative theories. 

As an application, we consider the evolution and effects 
of this fluid in cosmology. We show that there exists a 
class of scaling solutions which have an attractor solution 
with flxed equation of state Wfin- Off the attractor, this 
model possess an interesting dynamic where part of the 
energy density redshifts as dust while part of the energy 
density tracks any dominant background energy density. 
We use this curious property to construct a unifled dark 
energy/dark matter model where the limit = —1 
corresponds to standard ACDM. 



We also show that in this class of models, we can con- 
struct phantom models with wa^ < — 1 where there is 
no pathology when crossing the "phantom divide" at 
w — —1, at least classically. If we insist that the sys- 
tem satisfy the Null Energy Condition (i.e. w > —1), 
then we show that the kinetic term for perturbations is 
positive deflnite. 

Finally, we would like to conclude by emphasising that 
this class of theories provides a novel framework for cos- 
mological model building and exploring exotic states of 
matter. 
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Appendix A: Derivation of perturbed energy 
conservation 

The differential form of the perturbed constraint equa- 
tion (20) is 

6ip = (fil + ^i^i^^) 5ip + ^^, + 2<I>/^/i^ . (Al) 

Let us perturb the energy-conservation equation (9): 

(5u^V^e + 5e + 6e{e + p) + 3H {Se + dp) ^ . (A2) 

First of all, as usual we have 6u* — — (f> so that 

Suf'V^e = eSu'^ = . (A3) 

Further, we perturb the formula for the expansion Eq. 
(10) 

SO = fi^^s (Dip) - [fi'^nipfi^ + ^^^) Sip . 

For the perturbations of the d'Alambertian we have the 
standard result 

6 (aip) = -4^ip - 2$(^ - 6iJ$^ + 5tp + SHSip - ^Sip , 

which after using Eqs (11) and (20) can be written in 
terms of $ and Sep as 



S {Dip) = ^fil + fifi^^ + 3Hfi^ - Sip-3fi + H<i> 
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Therefore for the perturbation of the expansion we obtain 

Substituting this expression along with along the formula 
(A3) into Eq. (A2) results in 



which corresponds to the standard case (c.f. [36, p. 312, 
Eq. (7.105)]). 



[1] T. Padmanabhan and T. R. Choudhury, "Can the [15 
clustered dark matter and the smooth dark energy arise 
from the same scalar field?," Phys. Rev. D66 (2002) 
081301, arXiv:hep-th/0205055. [16 

[2] R. J. Scherrer, "Purely kinetic k-essence as unified dark 

matter," Phys. Rev. Lett. 93 (2004) 011301, [17 
arXiv : astro-ph/0402316. 

[3] D. Bertacca, S. Matarrese, and M. Pietroni, "Unified [18 
dark matter in scalar field cosmologies," Mod. Phys. 
Lett. A22 (2007) 2893-2907, arXiv: astro-ph/0703259. 

[4] A. Arbey, "Dark fiuid: A complex scalar field to unify 

dark energy and dark matter," Phys. Rev. D74 (2006) [19 
043516, arXiv:astro-ph/0601274. 

[5] A. Y. Kamenshchik, U. Moschella, and V. Pasquier, 
"An alternative to quintessence," Phys. Lett. B511 
(2001) 265-268, arXiv:gr-qc/0103004. [20 

[6] M. C. Bento, O. Bertolami, and A. A. Sen, "The revival 
of the unified dark energy - dark matter model?," Phys. 
Rev. D70 (2004) 083519, arXiv:astro-ph/0407239. [21 

[7] S. Capozziello, S. Nojiri, and S. D. Odintsov, "Unified 
phantom cosmology: infiation, dark energy and dark 
matter under the same standard," Phys. Lett. B632 [22 
(2006) 597-604, arXiv:hep-th/0507182. 

[8] S. Nojiri and S. D. Odintsov, "Unifying phantom 

infiation with late-time acceleration: Scalar [23 
phantom-non-phantom transition model and generalized 
holographic dark energy," Gen. Rel. Grav. 38 (2006) 
1285-1304, arXiv:hep-th/0506212. [24 

[9] C. Quercellini, M. Bruni, and A. Balbi, "Affine equation 
of state from quintessence and k-essence fields," Class. 
Quant. Grav. 24 (2007) 5413-5426, arXiv : 0706 . 3667 
[astro-ph] . [25 
[10] E. V. binder and R. J. Scherrer, "Aetherizing Lambda: 
Barotropic Fluids as Dark Energy," Phys. Rev. D80 

(2009) 023008, arXiv: 0811. 2797 [astro-ph]. [26 
[11] O. F. Piattella, D. Bertacca, M. Bruni, and 

D. Pietrobon, "Unified Dark Matter models with fast 
transition," JCAP 1001 (2010) 014, arXiv : 0911 . 2664 [27 
[astro-ph. CO] . 
[12] C. Gao, M. Kunz, A. R. Liddle, and D. Parkinson, 

"Unified dark energy and dark matter from a scalar [28 
field different from quintessence," Phys. Rev. D81 

(2010) 043520, arXiv: 0912. 0949 [astro-ph. CO] . 

[13] P. Creminelli, G. D'Amico, J. Norena, L. Senatore, and [29 
F. Vernizzi, "Spherical collapse in quintessence models 
with zero speed of sound," JCAP 1003 (2010) 027, [30 
arXiv: 0911. 2701 [astro-ph. CO] . 

[14] E. Komatsu et al, "Seven- Year Wilkinson Microwave 

Anisotropy Probe (WMAP) Observations: Cosmological [31 
Interpretation," arXiv: 1001 . 4538 [astro-ph. CO] . 



G. N. Felder, L. Kofman, and A. Starobinsky, "Caustics 

in tachyon matter and other Born-Infeld scalars," JHEP 

09 (2002) 026, arXiv:hep-th/0208019. 

A. Sen, "Tachyon matter," JHEP 07 (2002) 065, 

arXiv:hep-th/0203265. 

A. Sen, "Rolling Tachyon," JHEP 04 (2002) 048, 
arXi V : hep-th/02032 1 1 . 

N. Arkani-Hamed, H.-C. Cheng, M. A. Luty, 

S. Mukohyama, and T. Wiseman, "Dynamics of Gravity 

in a Higgs Phase," JHEP 01 (2007) 036, 

arXiv : hep-ph/0507120. 

N. Arkani-Hamed, H.-C. Cheng, M. A. Luty, and 
S. Mukohyama, "Ghost condensation and a consistent 
infrared modification of gravity," JHEP 05 (2004) 074, 
arXiv:hep-th/0312099. 

S. Mukohyama, "Caustic avoidance in Horava-Lifshitz 
gravity," JC^P 0909 (2009) 005, arXiv: 0906 . 5069 
[hep-th] . 

D. Bias, O. Pujolas, and S. Sibiryakov, "On the Extra 
Mode and Inconsistency of Horava Gravity," JHEP 10 
(2009) 029, arXiv: 0906. 3046 [hep-th]. 

P. Horava, "Quantum Gravity at a Lifshitz Point," 
Phys. Rev. D79 (2009) 084008, arXiv : 0901 . 3775 
[hep-th] . 

J. D. Brown and K. V. Kuchar, "Dust as a standard of 
space and time in canonical quantum gravity," Phys. 
Rev. D51 (1995) 5600-5629, arXiv:gr-qc/9409001. 
S. M. CarroU, M. Hoffman, and M. Trodden, "Can the 
dark energy equation-of-state parameter w be less than 
-1?," Phys. Rev. D68 (2003) 023509, 
arXiv : astro-ph/0301273. 

S. Dubovsky, T. Gregoire, A. Nicolis, and R. Rattazzi, 
"Null energy condition and superluminal propagation," 
JHEP 03 (2006) 025, arXiv:hep-th/0512260. 
R. P. Woodard, "Avoiding dark energy with 1/R 
modifications of gravity," Lect. Notes Phys. 720 (2007) 
403-433, arXiv:astro-ph/0601672. 
J. M. Cline, S. Jeon, and G. D. Moore, "The phantom 
menaced: Constraints on low-energy effective ghosts," 
Phys. Rev. D70 (2004) 043543, arXiv :hep-ph/0311312. 
R. Emparan and J. Garriga, "Non-perturbative 
materiahzation of ghosts," JHEP 03 (2006) 028, 
arXiv:hep-th/0512274. 

T. Jacobson and D. Mattingly, "Einstein-aether waves," 
Phys. Rev. D70 (2004) 024003, arXiv:gr-qc/0402005. 
S. M. Carroll and E. A. Lim, "Lorentz-violating vector 
fields slow the universe down," Phys. Rev. D70 (2004) 
123525, arXiv: hep-th/ 0407149. 

E. A. Lim, "Can we see Lorentz-violating vector fields 
in the CMB?," Phys. Rev. D71 (2005) 063504, 



13 



arXiv : astro-ph/0407437. 

[32] B. F. Schutz, "Perfect fluids in general relativity: 

Velocity potentials and a variational principle," Phys. 

Rev. D 2 (1970) no. 12, 2762-2773. 
[33] C. Armendariz-Picon, T. Damour, and V. F. 

Mukhanov, "k-inflation," Phys. Lett. B458 (1999) 

209-218, arXiv:hep-th/9904075. 
[34] C. Armendariz-Picon, V. F. Mukhanov, and P. J. 

Steinhardt, "Essentials of k-essence," Phys. Rev. D63 

(2001) 103510, arXiv:astro-ph/0006373. 
[35] C. Armendariz-Picon, V. F. Mukhanov, and P. J. 

Steinhardt, "A dynamical solution to the problem of a 

small cosmological constant and late-time cosmic 

acceleration," Phys. Rev. Lett. 85 (2000) 4438-4441, 

arXiv : astro-ph/0004134. 
[36] V. Mukhanov, Physical foundations of cosmology. 

Cambridge University Press, 2005. 
[37] J. Garriga and V. F. Mukhanov, "Perturbations in 

k-inflation," Phys. Lett. B458 (1999) 219-225, 

arXiv:hep-th/9904176. 
[38] R. R. CaldweU, "A Phantom Menace?," Phys. Lett. 

B545 (2002) 23-29, arXiv:astro-ph/9908168. 
[39] A. Vikman, "Can dark energy evolve to the phantom?," 

Phys. Rev. D71 (2005) 023515, 

arXiv : astro-ph/0407107. 
[40] R. R. Caldwell and M. Doran, "Dark-energy evolution 

across the cosmological-constant boundary," Phys. Rev. 

D72 (2005) 043527, arXiv:astro-ph/0501104. 
[41] W. Hu, "Crossing the phantom divide: Dark energy 

internal degrees of freedom," Phys. Rev. D71 (2005) 

047301, arXiv:astro-ph/0410680. 
[42] G.-B. Zhao, J.-Q. Xia, M. Li, B. Feng, and X. Zhang, 

"Perturbations of the Quintom Models of Dark Energy 

and the Effects on Observations," Phys. Rev. D72 

(2005) 123515, arXiv:astro-ph/0507482. 

[43] A. A. Sen, "Reconstructing K-essence," JCAP 0603 

(2006) 010, arXiv:astro-ph/0512406. 

[44] L. R. Abramo and N. Pinto-Neto, "On the Stability of 

Phantom K-essence Theories," Phys. Rev. D73 (2006) 

063522, arXiv:astro-ph/0511562. 
[45] M. Kunz and D. Sapone, "Crossing the Phantom 

Divide," Phys. Rev. D74 (2006) 123503, 

arXiv : astro-ph/0609040. 
[46] E. Babichev, V. Mukhanov, and A. Vikman, 

"k-Essence, superluminal propagation, causality and 

emergent geometry," JHEP 02 (2008) 101, 

arXi V : 0708 . 056 1 [hep-th] . 
[47] B. Feng, X.-L. Wang, and X.-M. Zhang, "Dark Energy 

Constraints from the Cosmic Age and Supernova," 



Phys. Lett. B607 (2005) 35-41, 

arXiv : astro-ph/0404224. 
[48] Y.-F. Cai, E. N. Saridakis, M. R. Setare, and J.-Q. Xia, 

"Quintom Cosmology: theoretical implications and 

observations," arXiv : 0909 . 2776 [hep-th]. 
[49] H. Zhang, "Crossing the phantom divide," 

arXiv: 0909. 3013 [astro-ph.CO] . 
[50] E. J. Copeland, M. Sami, and S. Tsujikawa, "Dynamics 

of dark energy," Lnt. J. Mod. Phys. D15 (2006) 

1753-1936, arXiv:hep-th/0603057. 
[51] B. Boisseau, G. Esposito-Farese, D. Polarski, and A. A. 

Starobinsky, "Reconstruction of a scalar-tensor theory 

of gravity in an accelerating universe," Phys. Rev. Lett. 

85 (2000) 2236, arXiv:gr-qc/0001066. 
[52] R. Gannouji, D. Polarski, A. Ranquet, and A. A. 

Starobinsky, "Scalar-tensor models of normal and 

phantom dark energy," JCAP 0609 (2006) 016, 

arXiv : astro-ph/0606287. 
[53] W. Hu and I. Sawicki, "Models of f(R) Cosmic 

Acceleration that Evade Solar-System Tests," Phys. 

Rev. D76 (2007) 064004, arXiv : 0705 . 1158 

[astro-ph] . 

[54] H. Motohashi, A. A. Starobinsky, and J. Yokoyama, 
"Phantom boundary crossing and anomalous growth 
index of fluctuations in viable f(R) models of cosmic 
acceleration," arXiv : 1002 . 1141 [astro-ph.CO]. 

[55] M. Park, S. Watson, and K. M. Zurek, "A Unified 
Approach to Cosmic Acceleration," arXiv : 1003 . 1722 
[hep-th] . 

[56] P. Creminelli, M. A. Luty, A. Nicolis, and L. Senatore, 

"Starting the universe: Stable violation of the null 

energy condition and non-standard cosmologies," JHEP 

12 (2006) 080, arXiv:hep-th/0606090. 
[57] P. Creminelli, G. D'Amico, J. Norena, and F. Vernizzi, 

"The Effective Theory of Quintessence: the vu < —1 

Side Unveiled," JCAP 0902 (2009) 018, 

arXiv: 0811. 0827 [astro-ph]. 
[58] C. Deffayet, O. Pujolas, I. Sawicki, and A. Vikman, 

"Imperfect dark energy from kinetic gravity braiding," 

in preparation (2010) . 
[59] S. Nojiri and S. D. Odintsov, "Inhomogeneous equation 

of state of the universe: Phantom era, future singularity 

and crossing the phantom barrier," Phys. Rev. D72 

(2005) 023003, arXiv:hep-th/0505215. 
[60] S. Capozziello, V. F. Cardone, E. Ehzalde, S. Nojiri, 

and S. D. Odintsov, "Observational constraints on dark 

energy with generalized equations of state," Phys. Rev. 

D73 (2006) 043512, arXiv:astro-ph/0508350. 



